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Direct study of the optical H vibrations of diluted Pd-H solid solutions by inelastic neutron 
scattering (INS) is impeded due to the very small limiting hydrogen concentration, x max < 0.001 
at low temperatures. We report study of disordered Pdo.8Auo.2Ha; solutions with concentrations 
varying from x — 0.03 to 0.74. Our INS investigation showed that the fundamental optical H 
band in both concentrated and dilute Pdo.8Auo.2Hj; solutions consists of a sharp peak with a broad 
shoulder towards higher energies. We argue that the optical band in infinitely diluted H solutions 
in Pd consists of a peak with a pronounced shoulder too. Implications for the isotopic effects in the 
hydrogen solubility in Pd are considered in the present paper. 



Introduction. In the fee lattice of palladium, 
hydrogen occupies octahedral interstitial positions 
with cubic symmetry. One could therefore expect 
that in diluted solid H solutions in Pd the funda- 
mental band of optical H vibrations should be re- 
duced to a narrow line of 3-fold degenerate, non- 
interacting isotropic local oscillators. A few at- 
tempts to test this assumption by inelastic neu- 
tron scattering (INS) gave ambiguous results due 
to the very small limiting hydrogen concentration, 
Xmax < 0.001, of the PdHz solutions at low temper- 
atures~0, 0] • 

We increased x max by alloying Pd with 20 at% Au 
and studied disordered Pdo.sAuo.2H1; solutions with 
concentrations varying from x = 0.03 to 0.74. The 
samples were prepared at hydrogen pressures from 
1 bar to 75 kbar. According to 197 Au Mossbauer 
studies 3, hydrogen atoms in the Pd0.sAu0.2Ha; 
solutions with x < (0.8) 6 « 0.26 can only occupy 
interstices having no Au neighbours. This abates 
the effect of H-Au interactions in the solutions with 
lower H concentrations and effectively prevents 
hydrogen clustering at low temperatures. 

Inelastic neutron scattering studies. Our 

INS investigation at 5 K (INl-BeF spectrometer, 
Institute Laue-Langevin, Grenoble) showed that the 
fundamental optical H band in both concentrated 
and diluted Pd0.sAu0.2Hz solutions consists of a 
sharp peak with a broad shoulder towards higher 
energies. Fig. [1] compares the optical spectra 
for the stoichiometric PdH sample 3] and the 
Pdo.8Auo.2Ho. 03 sample with the minimum H 
concentration. The shoulder amounted to approx- 
imately 0.44 integral intensity of the main peak in 
PdH and demonstrated no tendency to vanish with 
decreasing H concentration in the Pd0.sAu0.2Ha; 
samples. 

Both gold and hydrogen expand the lattice of pal- 
ladium. If the position Hloq of the main optical peak 
of the PdH and Pdo.sAuo.2H2; solutions is plotted 
as a function of the atomic volume, its linear ex- 
trapolation to the volume of pure Pd metal gives 
Hluq = 69.1 meV. This well agrees with an experi- 
mental value of 68.5 ± 0.5 meV for PdH .oo2 111 and 
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Fig. 1. Experimental one-phonon generalized 
vibrational densities of states, G(fvjS), corrected for 
the Debye- Waller factor. Open circles - fee PdH 
measured at 25 K Q. Solid circles - 
fee Pdo.8Auo.2Ho. 03 measured at 5 K with the 
INl-BeF spectrometer at ILL, Grenoble. 



therefore shows that the extrapolation procedure is 
meaningful. The extrapolation also shows that the 
spectrum Pd0.sAu0.2H0. 03 should well represent the 
fundamental optical band of the infinitely diluted 
Pd(H) solution. The position of the centre of grav- 
ity of this band gives the weight-averaged energy 
ftwaver ~ 80 meV for optical vibrations in Pd(H). 

One of the possible explanations of the high- 
energy shoulder in the spectra of diluted hydrogen 
solutions in palladium is the smearing of the peak 
due to what can be classically interpreted as the 
changes in the high-frequency vibrations of the H 
atom in a slowly varying "instant" potential of the 
slowly moving heavy metal atoms. Mention in this 
connection that a molecular-dynamics calculation 
for an H atom in the Pd matrix Q gave a Fourier 
transform of the time correlation function that qual- 
itatively reproduces the experimental peak with a 
broad shoulder in the INS spectra of the diluted Pd- 
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Fig. 2. Available experimental data on the 
solubility of hydrogen isotopes in Pd in the limit of 
infinite dilution. 

Au-H solutions. 



H2 gas and Gg r> (T) is the Gibbs energy of 1 mol of 
H atoms in Pd at infinite dilution. Subtracting this 
equation from an analogous equation for the Pd-D 
solutions yields the expression: 

AG DH /RT = [Gg(T) - G£?(T)] /RT = 
\ In [K (T)/Kk(T)] + \ [G D2 - Gh 2 ] /RT (1) 

The standard Gibbs energy for the gaseous phase 
of each hydrogen isotope has been accurately deter- 
mined and tabulated and is available in the litera- 
ture @. It is therefore advisable to calculate the 
right-hand side of equation (JTJ) and to consider it 
as the experimental AG'^/RT dependence to be 
fitted by the difference [Gg (T) - G£(T)} /RT that 
only depends on the properties of the Pd(D) and 
Pd(H) solutions. The temperature dependences of 
AG'^/RT and also AG^/RT are presented in 
Fig. [3] by symbols. 

Due to the large ratio of the atomic masses of 
hydrogen and palladium (mpd = 106.4 u), the dif- 
ference Gg' (T) — Gg 3 (T) is nearly solely determined 
by the difference in the free energies of optical vibra- 
tions of D and H atoms, the energies being counted 
from the bottom of the potential well for these atoms 
in the Pd(D) and Pd(H) solutions, correspondingly 
H and D atoms in Pd metal were earlier con- 
sidered as three-dimensional Einstein oscillators be- 
cause the main optical peaks had only been observed 
in the INS spectra of diluted PdR x and PdD^ solu- 
tions In the Einstein model: 



Implications for the isotopic effects in the hy- 
drogen solubility in palladium. The weight- 
averaged energy hui avcr of optical hydrogen vibra- 
tions plays an important role in the thermodynamics 
of the Pd-H system. The value Hloq = 68.5 meV of 
the main peak [l[ mistakenly used in the calculations 
instead of Sw aver « 80 meV since 1976 has caused 
much difficulty in the interpretation of experimental 
results and, in particular, led to the ill-founded con- 
clusion that the optical hydrogen vibrations in Pd 
are strongly anharmonic. 

Assuming that the spectrum of Pd0.sAu0.2H0. 03 
depicted in Fig. [1] represents the phonon density of 
states, g(ui), of the infinitely diluted Pd(H) solution, 
we reanalyzed the experimental results in literature 
concerning the protium, deuterium and tritium sol- 
ubility in palladium and demonstrated that these 
results agree with the harmonic model. 

The literature results are shown in Fig. [2] in 
the form of temperature dependences of the Siev- 
erts constant determined as — [H2] / [H] 2 = 
Pn 2 (l — x) 2 /x 2 for the investigated diluted PdH K 
solutions and in a similar way for the PdD^ and 
PdT^ solutions. In the case of infinitely diluted so- 
lutions, one can write: 

lnK H (T) = (2Gg>(T) - Gh 2 (T,P ))/RT, 

where Gn 2 (T, Pq) is the standard Gibbs energy of 



[G^(T)-G^(T)} /RTk 

+ ^(hw D -hiu H )/kT + A lm /kT (2) 

where Hlj is the energy of the first harmonics, k is 
the Boltzmann constant and Adh is the difference in 
the bottom energies of the potential wells for D and 
H atoms. The usage of huJn = 68.5 meV and fku^ = 
48 meV from ref. [1( gives the dependence shown in 
Fig. [3 by the thin solid line labeled as "D-H [1]". 
The line irrecoverably disagrees with experiment. 

If the densities of optical phonon states, <7h(w) and 
gi)(u), are known for the infinitely diluted Pd(H) 
and Pd(D) solutions, one can write 

\G<g{T)-G%{T)]/RTK 




ln[l — exp(— huj/kT)] \gn(uj) — g-n(uj)]dui+ 



3 f°° 

-J hw[gu{w) — gn(u>)]duj / kT + A^/kT, (3) 

where the g(uj) spectra are normalized so that 
/ g(u))dus = 1. 

Assuming that gn{^) is proportional to the gen- 
eralized density of states for the Pd0.sAu0.2H0. 03 
sample shown in Fig. [T]and that <?d(w) = rgu_(ruj), 
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optical vibrations in Pd(H). 
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Fig. 3. The dimensionless differences AGdh/RT 
and AGtd /RT calculated using the experimental 

data from Fig. [5] and results of fitting these 
differences by modelling the dynamics of hydrogen 
isotopes in Pd metal (solid and dashed lines, see 
text). 

we fitted AG^ / RT with equation (J3J) using r and 
Adh as fitting parameters. The resulting depen- 
dence is shown in Fig. [3] by the thick solid line la- 
beled "D-H for g(uj)" . As one can see, it well agrees 
with experiment. The obtained optimum value of 
r = V2.06 ~ yj tod / wh suggests the nearly har- 
monic character of optical vibrations in Pd(H) and 
Pd(D). The optimum value of Adh = —3.3 meV 
indicates that the bottom of the potential well for 
H atoms in Pd lies deeper than for D atoms by 
3.3 meV. 

Fitting the AG^/RT from ref. Q in a similar 
way yields a rather good approximation to exper- 
iment (upper thick solid line in Fig. [3]) with r = 
VF52 w y'mT/mD and A TD ~ 0. The AG C ^/RT 
dependence from ref. [|[ (open stars in Fig. [3J is less 
accurate than that from ref. |7| and it cannot be fit- 
ted with physically acceptable values of the fitting 
parameters. 

Interestingly the AG D X ^ / RT and AG^/RT de- 
pendences can be well approximated with equation 
((2} (see dashed lines in Fig. [3]) using virtually the 
same values of r and A as above, if hwu is set equal 
to the weight-averaged energy fejj avC r = 80 meV of 



Conclusions. The inelastic neutron scattering in- 
vestigation of the Pd-Au-H solutions has for the first 
time given the opportunity to establish the spectrum 
of optical vibrations in diluted Pd-H solutions. 

In contrast to the commonly held view, the band 
of optical H vibrations in the Pd-H solutions does 
not shrink to a narrow line of 3-fold degenerate local 
oscillators on infinite dilution. Instead, the funda- 
mental optical peak has a broad and intense shoul- 
der towards higher energies. As a result, the weight- 
averaged energy fiw a vcr of optical vibrations in di- 
luted Pd-H solutions is equal to about 80 meV and 
not to 68.5 meV as it was assumed earlier. 

The value of 7uJ a ver = 80 meV well describes all 
available experimental data on the H and D solu- 
bility in Pd in the harmonic approximation. This 
affords no basis for the widely spread opinion that 
optical vibrations in the diluted Pd-H solutions are 
strongly anharmonic. 

The analysis of the solubility of hydrogen isotopes 
in palladium shows that the bottom of the potential 
well is located at approximately the same energy 
for the D and T atoms and lies by some 3.3 meV 
deeper for the H atoms. The force constants are a 
little higher for the H atoms than for the D atoms 
and coincide for the D and T atoms within the 
experimental error. 
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